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Abstract 

O/ Efficient link scheduling in a wireless network is challenging. Typical optimal algorithms require 

solving an NP-hard sub-problem. To meet the challenge, one stream of research focuses on finding 
^ ' simpler sub-optimal algorithms that have low complexity but high efficiency in practice. In this paper, 

G^ ■ 

Q^ , we study the performance guarantee of one such scheduling algorithm, the Longest-Queue-First (LQF) 

^H ■ 

^f^ • algorithm. It is known that the LQF algorithm achieves the full capacity region. A, when the interference 

C^^ ' graph satisfies the so-called local pooling condition. For a general graph G, LQF achieves (i.e., stabilizes) 

o : 

a part of the capacity region, a*{G)A, where (7*{G) is the overall local pooling factor of the interference 

graph G and (t*{G) < 1. It has been shown later that LQF achieves a larger rate region, I]*(G)A, 

. ■ where I]*(G) is a diagonal matrix. The contribution of this paper is to describe three new achievable 

X ■ 

^ ' rate regions, which are larger than the previously-known regions. In particular, the new regions include 



all the extreme points of the capacity region and are not convex in general. We also discover a counter- 
intuitive phenomenon in which increasing the arrival rate may sometime help to stabilize the network. 
This phenomenon can be well explained using the theory developed in the paper 

Index Terms 
Wireless Networks Scheduling, Longest Queue First PoUcy, Stability, Local Pooling, Interference 



I. Introduction 

One of the long-standing challenges for wireless networks is how to utilize the communication medium 
efficiently when links interfere with each other. This paper is primarily concerned with an interference 
model called the protocol model, where two wireless links that interfere with each other are prohibited 
to transmit data simultaneously |[T1^ 121. El- For the protocol model, a scheduling algorithm strives to 
select a set of non-interfering links for transmission in every time slot. 

Finding efficient schedules can be very difficult. Tassiulas and Ephremides lH showed that if the queue 
sizes for the links (which are nodes in the interference graph) are viewed as weights and a maximum 
weight independent set (MWIS) of the interference graph is selected as the schedule in each time slot, then 
the queues of the wireless network can be stabilized for any arrival rate vector inside the capacity region. 
However, finding an MWIS is NP-hard in general. Even under the more restricted k-hop interference 
model, finding an MWIS is still NP-hard for A; > 2 Q, IS. For the 1-hop interference model, the problem 
of finding an MWIS reduces to maximum weight matching and the complexity is 0(|y|^), where \V\ is 
the number of wireless links IS. Hence, scheduling using MWIS is inapplicable to large networks. 

To reduce the complexity, some simple sub-optimal scheduling algorithms have been introduced Q, 
11 , HI, Uni, 0, ini, ini- The Longest Queue First (LQF, also known as the greedy maximal 
schedule) policy is recognized for its high performance in practice ||6l. The LQF schedule chooses links 
in a decreasing order of the queue sizes while conforming to the interference constraints. Dimakis and 
Wakand showed that the LQF algorithm achieves (stabilizes) the entire interior of the capacity region, 
A, when the interference graph G satisfies the so-called overall local pooling condition fT4l . For general 
cases, Joo et al. introduced a parameter called the overall local pooling factor, denoted by a*{G), where 
< (J*{G) < 1, based on the topology of the interference graph G; they showed the LQF algorithm 
achieves a subset of the capacity region, a*{G)A 13. Several other authors studied how to check the 
local pooling condition or estimate cr*{G) for specific graphs ITSl . lT6ll . ITTl . 

A single -parameter performance characterization of LQF suggests a uniform rate reduction on the 
links. However, it is possible that the links are subject to heterogeneous interference relations and some 
links can perform better than others. To capture the performance heterogeneity, a multiple-parameter 
characterization of the stabilizable region by LQF was established in ITSl . It was shown that LQF can 
achieve a larger rate region, S*(G)A, where i;*(G) is a diagonal matrix. Each diagonal entry of S*(G) 
corresponds to a link and it summarizes the link's interference constraints. 

Even this multiple-parameter characterization of LQF underestimates the stability region. For instance. 



it excludes some parts of the capacity region that are obviously stabilizable by LQF. To progress further 
toward complete performance characterization, there is a need to go beyond the current framework of 
linear transformations on the capacity region. The goal of this study is to establish such a "non-linear" 
framework and expand our knowledge about the achievable rate region by LQF. The main contribution is 
to describe three new achievable rate regions (il, Ac and Ar), which are all larger than the previously- 
known regions. More precisely, we show that S*(G)A C Q, and 0,° C Ac ^ A/j. Furthermore, the 
closures of the new regions include all the extreme points of the capacity region and are not convex 
in general. This is in contrast to previously-known regions of stability, which are all convex and, in 
general, exclude some extreme points of the capacity region because they each are derived by reducing 
the capacity region through a linear transformation. We show that the new regions of stability (or their 
closures) are convex if and only if they are identical to the capacity region itself. The result implies 
that, when LQF cannot achieve the full capacity region, the largest achievable region, which is yet to be 
discovered, cannot be convex. 

The characterization of the LQF performance has been substantially improved with these new stability 
regions. For instance, we have found that, for an arbitrarily large A; > 0, there are cases where an arrival 
rate vector A is outside all the previously-known stability regions but kX is in Q,. In other words, the 
previously-known stability regions can underestimate the performance of LQF by an arbitrarily large 
factor in certain cases, whereas the new regions can avoid such poor estimates. 

The study has also yielded an interesting, counter-intuitive finding that increasing the arrival rates may 
sometime help to stabilize the network. We have discovered an example where a rate vector achievable 
by LQF point-wise dominates another rate vector not achievable by LQF. It turns out the former vector 
is in the stability region Ac whereas the latter is not. 

We next summarize the key ideas of the paper. Our theory is developed based on considering the fluid 
limit of an unstable network. A typical scenario is that the maximum queue size has an overall trend to 
grow indefinitely, which requires that, at some time t, a subset of the current longest queues continues 
to grow. From the set of the longest queues at time t, there is a subset that grows at the fastest rate 
and remains the longest in the next infinitesimal time interval. Denote this subset by S. Under LQF, the 
queues in S will be served with priority in the next small time interval, which implies that the average 
service rate vector, when restricted to S, comes from the convex hull of the maximal schedules with 
respect to S. This convex hull is denoted by Co{Ms)- For the queues in S, the arrival rates must be 
larger than the service rates. The discussion motivates the definition of a strictly dominating vector for a 
queue set S, which is a vector A, when restricted to S, strictly dominating at least one vector in Co{Ms)- 



After removing the union of the strictly dominating vectors, where the union is over all possible subsets 
of the queues, we get $7. 

Key to the development about Ac is a refinement to the notion of strictly dominating vectors, which 
is called uniformly dominating vectors. For the aforementioned queue set S, the arrival rates not only 
must be larger than the service rates, but also larger by the same amount, so that the queues in S grow 
at the same rate. The removal of all the uniformly dominating vectors gives Ac. Although the closure of 
Ac contains Q,, there is value in studying and reporting the results about both regions. First, the theory 
about Q provides building blocks for proving some of the results about Ac. Second, appears to be 
well connected to the notion of local poohng in |[T4l . thus, providing some continuity in the theoretical 
development, whereas Ac does not appear so. 

Throughout, we assume i.i.d. and mutually independent arrival processes. As Dimakis and Walrand 
pointed out, for the same average arrival rate vector, whether the arrival processes have zero or non-zero 
variances leads to significantly different stability behavior (the former is the case of deterministic arrivals 
with constant rates) |[T4l . They established a queue separation result for the case of non-zero variances 
and developed a rank condition that leads to queue separation. We generalize the rank condition. Then, 
we extend Ac to a larger stability region Aji for the case of non-zero variances. We also show the 
closures of Ac and A/j are the same. 

Finally, we relate the problems of finding stability conditions under LQF to several problems in the 
fractional graph theory |[T9l . The latter provide tools for studying the stability regions introduced by the 
paper and for characterizing the set a-local pooling factor given in |[T8l . 

The rest of the paper is organized as follows. In Section |lll we specify the models and notations. 
In Sections |lll] and |lVl we introduce the il and A (Ac and A/j) regions, respectively. In Section |Vl 
we introduce the fractional coloring and related problems that are relevant to the study of the stability 
regions. In Section |Vll we give some simulation results to confirm aspects of the theory. The conclusion 
is given in Section |VII[ 

II. Preliminaries 

In our model, a wireless network is represented by an undirected interference (or conflict) graph 
G = {V,E), where the node set V represents the set of physical, wireless links in the network and the 
edge set E represents the interference relation among the physical links. Two nodes in G are connected 
with an edge whenever the physical links they represent interfere with each otherj^ 

'All the graphs in this paper are interference graphs, unless specified otherwise. 






We assume the node 



set V is arbitrarily indexed from 1 to \V\, and hence, V can be written asF = {l,...,|F|}. 

Given a subset of nodes 5 C y, we denote Gs = [S, L) to be the subgraph of G induced by the 
nodes in S. In other words, an edge (u, v) belongs to L if and only if u, w G S" and (u, v) G E. 

We assume a time-slotted system. The capacity of each wireless link is normalized to 1 per time slot. 
There is a queue associated with each wireless link at the transmitter. We assume single-hop traffic. Traffic 
arrives at the transmitter side of a link, joining the queue and waiting for transmission; after transmission, 
it leaves the network. We assume i.i.d. and mutually independent arrival processes to the queues. It is 
easy to see that, under the LQF schedule with either deterministic or typical random tie-breaking rules, 
the joint queue process is Markovian. By stability, we mean the Markov process is positive recurreno 

A schedule is denoted by a |y [-dimensional 0-1 vector, where a value 1 in an entry indicates the 
corresponding link is active and otherwise. A schedule is feasible if and only if the links that are 
active do not interfere with each other. A feasible schedule is said to be maximal if no additional links 
can be activated without violating the interference constraints. Therefore, every feasible schedule is an 
independent set of G and every maximal schedule is a maximal independent set of G. 

For the graph G = {V,E), let My denote the set of all the maximal schedules and let Go{Mv) 
denote the convex hull of all the maximal schedules. When relevant, we also consider My to be the 
matrix whose columns are all the maximal schedules, with arbitrary indexing of the schedules. Similarly, 
for a node-induced subgraph Gs = {S,L), let Ms be the set (or matrix) representing all the maximal 
schedules of Gs and let Co{Ms) be the convex hull of all the maximal schedules in Ms- 

The capacity region A of a network is defined as the set of arrival rate vectors that are supportable by 
time sharing of the feasible schedules. Equivalently, 

A = {A I < A < /i for some /x G Go{My)]. (1) 

For a non-empty subset of nodes S* C F, the capacity region is defined analogously by replacing Go{My) 
with Go{Ms) in ([T|) and is denoted by A5. In the above, X < ji means that vector A is component-wise 
less than or equal to vector /i. The interior of the capacity region can be written as 

A° ={A |0 < \<^l for some /i G Co{My)}. (2) 

The interior of the capacity region thus defined can be stabilized by the MWIS schedule and any rate 
vector outside the capacity region cannot be stabilized by any schedule H. 

^Without loss of generality, we assume the Markov Chain is irreducible. See (4) for general cases. 



Given a |F| -dimensional vector A, the |5| -dimensional vector [X]s represents the restriction of A to 
the set S C V. That is, [A] 5 contains only those components of A which correspond to the nodes in 5. 

For a vector fi defined for a node set, let fii or ^{l) denote the component associated with / G V. Note 
that, if ;U G M.[^ , then the notation indicates the Ith component of fi. However, if ^u G m!^ for some 
non-empty S (^ V, then for I ^ S , fii or fi{l) is not necessarily the Ith component of /i. If jj. is any other 
type of vector, ^j denotes the ith component of jj.. We use e to represent the vector (1,1,...,!)'. The 
dimension of the vector e depends on the context. 

The capacity region for the whole graph and the capacity region for the subset 5 C 1/ has the following 
relationship. 

Lemma 1: An arrival rate vector A G A if and only if for all non-empty S CV, [X\s G A5. Likewise, 
A G A° if and only if for all non-empty S C y, [A]^ G A^. 

Proof: Suppose A G A. Then, A < ^u for some fi G Co{Mv)- It is easy to see that, for any non- 
empty subset 5 C y, there must exist a vector v G Co{Ms) such that [ij]s < v. Then, [A]s < v. Hence, 
[X]s G A5. The other direction is true by taking S = V. The last statement of the lemma can be proved 
similarly. ■ 

Throughout, in the statements about rate regions that involve topological concepts such as open/close 
sets and the interior of a set, the space is assumed to be the set of non-negative real vectors, i.e., M._^_ . 
Also, in the set-complement operation for any rate region, the whole set is understood to be the non- 
negative real vectors. For a set y C W_^_ , we let Y° and V^ denote the interior and complement of Y, 
respectively. 

In the LQF schedule, the links with longer queues are activated at a higher priority than those with 
shorter queues, subject to the interference constraints. The following may be considered as a reference 
implementation of this schedule. First, one of the links with the longest queue is selected to be in the 
schedule; ties are broken with either an arbitrary deterministic rule or randomly. All links with which the 
selected link interferes are removed from further consideration. Then, the same selection process repeats 
over the remaining links yet to be considered until no links remain to be considered. 
Remark: The following is the key mathematical property about LQF that is used throughout. Suppose, 
at time t, a non-empty set 5 C y dominates V — S in the sense that, for any i G S and any j ^ V — S, 
the queue size of i is greater than that of j. Then, the schedule used at t must be maximal when restricted 
to S (i.e., with respect to G5). 



III. Stability Region Q. under LQF 

In this section, we introduce a notion of strictly dominating vectors and construct a region denoted by 
i7 based on this notion. The Jl region is larger than a*{G)A and S*(G)A, which have previously been 
shown to be stabilizable by the LQF policy. Unlike those previously-discovered regions of stability, the 
ri region includes all the extreme points of A and it is not convex in general. 

A. Review of Set, Link and Overall a-local Pooling 

Set cr-local pooling has been studied in ifTSl . It has many interesting properties and is related to (overall) 
cr-local pooling defined in O. 

Definition 1: Given a non-empty set of nodes S CV, the set a-local pooling factor for S, denoted 
by cjj, is given by 

a*g =sup{(T I afi ':f> z^, for all fi,iy € Co{Ms) } (3) 

= inf{(T I ayi > z^, for some /i,z^ G Co{Ms) }■ (4) 

It has been shown that the set a-local pooling factor is equal to the optimal value of the following 
problem. 

cj% = miner, subject to (Jfi> v, IJ.,^ £ Co{Ms)- (5) 

The link u-local pooling factor is defined as follows. 

Definition 2: The local pooling factor of a link I G V, denoted by a^ , is given by 

a* =sup{o-|o-^ ;2i V for all 5 C y such that / G 5, and all ^t, v G Co{Ms)} (6) 

= 'm.{{a\a^ > u for some S C y such that / G S, and some ^, z^ G Co{Ms)}- (7) 

Comparing the definitions of cr|. and a^, we have 

(T;* = min (j%. (8) 

{SQV I l&S} 

The overall a-local pooling factor of the graph G = {V, E) is 

a*{G) = mino"/*. 
lev 

Let the diagonal matrix I1*(G) be defined by S*(G) = diag (cr;*)/gv'- It has been shown that cr*(G)A 

and S*(G)A are both regions of stability under LQF ||6l ifTSl . with the latter containing the former. 



B. Strictly Dominating Vectors and Vi Region 

We first discuss some intuition that leads to the construction of the Q. region. When the network is 
unstable, a typical situation is that the size of the longest queues has an overall trend of increase, if one 
ignores the short-time fluctuations. This would not have occurred if, for any subset 5 C y, the arrival 
rate is strictly less than the service rate at each node in S. Here, we imagine S is the set of nodes with 
the longest queues for an extended period of time. Then, over that period of time, the schedule on each 
time slot must be maximal when restricted to S and, by time sharing of such maximal schedules, the 
(average) service rate must be in Co{Ms)- The discussion motivates us to define the notion of strictly 
dominating vectors for a subset of the nodes. 

Definition 3: Given a non-empty node set 5" C y, a vector A G M^ is a strictly dominating vector 
of S if [X]s > V for some v G Co{pIs)- The region composed with all the strictly dominating vectors 
of S is called the strictly dominating region of S and is denoted by J\.s- That is, 

n5 = {A G Mf ' I \\\s > ly, for some u G Co{Ms)}. 

For convenience, if 5 = 0, we assume Us = 0. 
We are often interested in the complement of ITS': 

n| = {A G M^^' I [\]s 1> y, for all v G Co{Ms)] 

= {A G m!^' I for every u G Co{Ms), there exists / G S such that \{l) < u{l)}. 

Definition 4: The i7 region is defined by 

n= f] Ws. 

scv 
Remark. A vector A is outside i7 if and only if A G 115 for some non-empty node set S. Also, when 

restricted to the components corresponding to the nodes in 5, 115 is an open set (it is a union of open 

sets). Hence, 17 is a closed set. It can also be helpful to think = (IJ5cy ^S')'^- 

Lemma 2: Suppose an arrival rate vector A satisfies A G U°. Then, for any non-empty subset S <^ V 

and any u G Co{Ms), there exists / G S such that A(/) + e^, < z/(/), where eo > is a constant 

independent of S, u and I. 

Proof: Since A G 0,°, we have A + ee G il for some small enough e > 0. Suppose the conclusion 

of the lemma is not true. That is, suppose for any e > 0, there exists a non-empty subset 5 C y and 

u G Co{Ms) such that A(/) + e > v{l) for all I G S. We can choose e satisfying < e < e. Then, 

[A + ee]5 > [A + ee]s > f- Hence, A + ee G lis, which implies A + ee fi by Definition |4l leading to a 

contradiction. ■ 



C. Performance Guarantee of LQF in Q Region 

Therorem 3: If an arrival rate vector A satisfies A G fl°, then, the network is stable under the LQF 
policy. 

The full proof requires replicating most of the arguments in lfT4l . In the following, we only highlight 
the part of the argument that needs modification. 

Sketch of Proof: Consider the fluid limit of the queue processes, denoted by {qi{t)}t>o for all I € V. 
For a fixed (and regular) time instance t, let S be the set of those longest queues whose time derivatives 
at t, qi{t), are the largest under a given LQF policy instance. The queues in S will remain the longest 
with identical length in the next infinitesimally small time interval. 

The service rate vector, when restricted to S, must belong to the set Co{Ms). Roughly, this is because 
S contains all the queues that are the longest and remain the longest in the near future, and hence, as 
remarked earlier, every LQF schedule being used must be a maximal schedule when restricted to S. 

Now imagine u is the service rate vector for the nodes in S at time t. Since A € 0°, by Lemma |2l 
there exists a link I G S such that A(/) + eo < iy{l) for some constant eo > 0. Then, ij{l) — X{1) > eo- 
Hence, at any time instance, each of the longest queues decreases at a positive rate no less than Eq. This 
is sufficient to conclude that the original queueing process is a positive recurrent Markov process (see 
1201), which means the queues are stable. ■ 

Next, we show Q, contains the previously-known regions of stability for LQF. 

Lemma 4: The following holds: S*(G)A C n. 

Proof: Consider any vector A G S*(G)A. Let 5 C y be an arbitrary non-empty node set. Let 
n G argmaxfcg5cr^„. Since A C S*(G)A, by Lemma [Hand (H), [X\s G cr^^s ^ o-^As- Then, there exists 
a vector /i G A5 such that [AJ^ < <Jgfi. According to Definition[T] cjj/i '^ v for any v G Co{Ms)- Hence, 
[AJs' ;z^ V for any v G Co{J\Is\ implying A G 11^. Since S is chosen arbitrarily, A G n5c\/n5 = $7 by 
Definition IH ■ 

D. Shape of Region 

The previously-known regions of stability under LQF, such as S*(G)A, are derived by reducing the 
capacity region through a linear transformation. Since the capacity region A is convex, each of these 
derived stability regions is also convex. In contrast, we will show that the shape of the Q region is not 
convex in general. Furthermore, when the previously-known regions are not identical to A, they exclude 
many, if not most, of the extreme points of A. We will show that [7 contains all the extreme points of A. 
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Lemma 5: The set of all the independent sets of the interference graph G, i.e., the set of all the feasible 
schedules, has a bijection to the set of all the extreme points of A. 

Note that we consider the empty schedule where no link is activated a trivial independent set. Lemma |5] 
establishes a connection between the graph topology and the geometry of A in a vector space. The proof 
can be found in the Appendix. 

Lemma 6: Suppose A is a vector corresponding to an independent set of the interference graph G. 
Then, X £ Q. 

Proof: [X]s is an independent set of the node-induced subgraph Gs for any non-empty S* C T/. 
Then, [A]5 '^ v for any v G Go(Ms)^ which implies A n^. Since S is arbitrary, we must have A G $7. 

■ 
Corollary 7: All the extreme points of the capacity region A belong to O. 

Proof: This is a result of Lemma |5] and Lemma [6] ■ 




Fig. 1. The six-cycle graph, C&. 

As an example, let G be the six cycle graph in Fig. [U The arrival rate vector A = (1,0, 1,0, 1,0)' 
corresponds to an independent set, and hence, A € il. However, we know that S*(G) = diag(2/3, 
2/3, 2/3, 2/3, 2/3, 2/3). As a result, A S*(G)A. The example shows that Q. can be strictly larger than 
S*(G)A. In the example, — S*(G)A contains not only the extreme points. For instance, one can check 
that, for A = (7/10, 1/10, 7/10, 1/10, 7/10, 1/10)', A G SI but A ^ S*(G)A. 

The next example shows that, the previously-discovered stability regions (t*{G)K and S*(G)A can 
underestimate the performance of LQF by an arbitrarily large factor in certain directions and in certain 
cases, whereas can avoid such poor estimates. 

Lemma 8: For any A; > 0, there exists an interference graph G = {V, E) and an arrival rate vector A 
such that A ^ S*(G)A, but k\ G f7. 

Proof: Consider the bipartite graph in Fig. |2] with N pairs of nodes, where A^ = 4 in this particular 
case. It is almost a complete bipartite graph except that every corresponding pair of nodes (such as nodes 
1 and 2) does not have an edge between them. It is easy to check that S*(G) = diag(2/Af, 2/A^, ..., 2/A^). 
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Therefore, the rate vector A = (2/iV + e,0,0, ...,0)', where e > 0, is not in S*(G)A. For any A; > 0, 
we can find a large enough A^ and a small enough e such that k{2/N + e) < 1. Then, we have kX = 

{k{2/N + e),0,0,...,Oy inn. ■ 



Fig. 2. A bipartite graph. 





Fig. 3. The fl region and other relevant regions. The largest convex polytope is A. The entire shaded region is Q, which is 
not a convex set. 



Though we cannot draw various regions in a high-dimensional vector space, it may still be helpful 
to make a highly simplified illustration with Fig. [3] The whole capacity region A is convex. The region 
S*(G)A is derived by scaling down the capacity region A using the diagonal matrix S*(G). This sort 
of scaling usually cuts off many or most extreme points of A. The newly defined stability region is 
a superset of S*(G)A and U contains all the extreme points of A. The figure makes the point that Q, is 
not convex in general. We next show U is convex if and only if it is equal to A. 

Lemma 9: The following statements are equivalent. 

1. f] is a convex. 

2. G is an overall local pooling graph. 

3. f] = A. 

Proof: First, we prove that statement 1 implies statement 2. Suppose G is not overall local pooling. 
We claim that there must exist a non-empty set S C y and ^, z^ € Co{Ms) such that jj. > u. Since G 
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is not overall local pooling, there exists a non-empty set S" C y such that cr^ < 1, which implies that 
there exist fi, u G Co{Ms) and a^fi > u, according to ^. li u > 0, we have the required set S, and 
/i, z^ G Co{Ms) with ^f > u. If not, let H = {I £ S\iy{l) > 0}. Because u £ Co{Ms) and [v\s^h = 0, 
it is easy to show [v]h G Co{Mh]^ Because f^L G Co{Ms) and H (1 S, there must exist /i G Co{Mh) 
such that /i > [/i]^. Then, a*gjl > <Jg[fi]H > ['^l// > 0. Thus, fl > [uJh and /i, [zv]^ G Co{Mh). By 
renaming if to be S, /i to be /x and [i/]j|/ to be u, we have the required set S and /U, z/ G Co{Ms) with 
/i > i^. 

Let A G M^ be an extended vector from ^u such that [X\s = f^ and [AJy^s' = 0. According to Definition 
[3]and|4l A fi. Since ^u G Co(M5), we can write /i = YJi=i Oiim^, where '^iOti = 1 and Oj > for 
all i, and ?7i* for i = 1, . . . , IM5I are all the maximal schedules with respect to S. For each i, let m* be 
a |F|-dimensional vector extended from m*, such that rh}{j) = m^{j) when j € S and rh}{j) = when 
J 5. Clearly, each m* corresponds to an independent set of G. Hence, by Lemma [6l m* G O for all i. 
Since A = X]l=r ct*"^* ^^id A i7, we conclude that Q. is not convex. 

Next, we show that statement 2 implies statement 3. Since G is an overall local pooling graph, 5]*(G) = 
/ (the identity matrix). By Lemma IH A C Jl. Hence, 17 = A. 

Finally, statement 3 implies statement 1 since A is convex. ■ 

Remark. Suppose, for a given interference graph, the LQF algorithm does not achieve the full interior 
of the capacity region. Lemma |9] implies that is not convex. Furthermore, since the closure of the full 
stability region of LQF (which is unknown) contains Q., it contains all the extreme points of the capacity 
region A. Hence, the closure of the full stability region of LQF cannot be convex either, and it cannot 
be characterized by any linear transformation of the capacity region. 

IV. Stability Region A under LQF 

In this section, we develop a notion termed as uniformly dominating vectors. It leads to a stability 
region Ac, which is a superset of il°. When the arrival processes are not constant, i.e., when the variances 
of the i.i.d. arrival processes are non-zero, we obtain a stability region A/j, which contains Ac. 

^Suppose we write Ms = {m'')\^f , where each m' is a maximal schedule with respect to S. We can represent v as 
V = Yli^i oi-ii'^'^ < where '^^ on = 1 and a^ > for all i. Since [v\s-h = 0, we have [m'^\s-H ~ for each i. It is clear 
that [Tn^]H corresponds to an independent set of Gh, the subgraph of G induced by H. Moreover, by the maximality of m' 
with respect to 5", if nC{j) = for some j G H, it must be that nrC{k) = 1 for some k £ H and j and k interfere with 
each other, i.e., {j, k) G E. Therefore, [m^]H must be maximal with respect to H. Hence, by [i/Jh = X^'=f ai[m^]H, we get 
Hh G Co{Mh). 
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A. Motivating Examples 

Example 1: We will first give an example to show that an arrival rate vector \ ^ Q. can sometime be 
stabilized by LQF. Hence, there is a region larger than Q. that captures the performance of LQF more 
precisely. The example also contains hints about how such a region can be defined. 

Consider the six cycle graph G in Fig.[T] There are exactly five maximal schedules: s^ = (1, 0, 1, 0, 1, 0)', 
s^ = (0,1,0,1,0,1)', s^ = (1,0,0,1,0,0)', s^ = (0,1,0,0,1,0)', s^ = (0,0,1,0,0,1)'. Suppose the 
arrival rate vector is A = (5/12 + e, 1/3 + e, 1/3 + e, 1/3 + e, 1/3 + e, 1/3 + e)', where e > is some small 

|e. Then, one can check that // = \s^ + \i 



enough constant. Let e = (1, 1, 1, 1, 1, 1)', ^ = \e and v = le. Then, one can check that fj, = Is^ + is^ 



and I' = ^s^ + ^s^ + |s^, which implies that /x, i^ G Co{Mv)- For < e < 1/12, z/ < A < /x. Hence, 
A G A° and A by Definition [3] and H 

Consider the fluid limit of the queue processes under LQF, denoted by {qi{t)}t>o^ for all / € V. 
For a fixed (regular) time instance t, let S be the set of those longest queues whose time derivatives 
at t, qi{t), are the largest. The queues in S will remain the longest with identical length in the next 
infinitesimally small time interval. Since A € A°, [X]s G A^ by Lemma [T] If S y^V, it is a fact that the 
node-induced subgraph Gs satisfies the local pooling condition ITSl . An argument similar to that in the 
proof of Theorem [3] shows that the queues in S all have a negative drift. 

The case of 5 = y is more subtle. Since only the maximal schedules of G are used during the 
aforementioned infinitesimally small time interval, we can assume that the service rate vector is 7 = 
X]j=i ais\ where X]j=i «« = 1 and Oj > for all i. In the fluid limit, qi{t) = Xi — ^i for I £ V. By 
assumption, qi{t) should be identical for all nodes / G V. However, one can check that it is impossible 
to find such 7 for the given A. Therefore, the case oi S = V would not have occurred, and only the 
case oi S y^ V needs to be considered. Hence, G is stable under LQF for the given A, according to the 
discussion for the S ^ V case. 

Example 2: Let A^ = 0.7(1/2 - e, 1/2 - e, 1/2 - e, 1/2 - e, 1/2 - e, 1/2 - e)' and A^ = (1/2 - e, 1/2 - 
2e, 1/2 - 2e, 1/2 - 2e, 1/2 - 2e, 1/2 - 2e)' and e = 10^^ goth A^ and A^ are outside O. Interestingly, 
although A^ < A^, A^ cannot be stabilized by LQF while A^ can. This has been verified by simulation 
experiments under constant arrivals. We will next develop a theory that provides a larger stability region 
and also can explain this counter-intuitive example. 



B. Uniformly Dominating Vector and Ac Region 

Definition 5: Given a non-empty node set S 
dominating vector of S if [\]s = u + de for some u G Co{Ms) and scaler d > 0. The region composed 



Definition 5: Given a non-empty node set 5 C V^, a vector A G mL is said to be a uniformly 
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with all the uniformly dominating vectors of S is called the uniformly dominating region of S and is 
denoted by r^. That is, 

r^ = {A G Mi|_ I [X]s = 1^ + de, for some v G Co{Ms) and some scalar d > 0}. 

By convention, if 5 = 0, we assume Ts = 0- 
Definition 6: The Ac region is defined by 

Ac = n ^s- 

scv 
Remark. Note that a vector A is outside Ac if and only if A € Ts' for some non-empty node set 5. 

Lemma 10: For any non-empty 5" C y, r5 is closed. Hence, Ac is open. 

Proof: Let B = {de\d > 0}, where e is IS*] dimensional, and let C = Co{Ms). It is easy to see 
C is compact and B is closed. From Definition [51 Ts is B + C extended to the |y [-dimensional space. 
It can be shown that B + C is closed, and hence, Ts is closed. Then, Ac = flscy ^s ^^ open (with 
respect to the metric space M\^ ). ■ 

Lemma 11: Suppose A € Ac and suppose 5 C 1/ is a non-empty node set. If u — [X]s = de for some 
u G Co{Ms), then d > eo, for some eo > independent of S and u. 

Proof: Suppose v — [X\s = de (here, e is of jSI -dimension) for some v G Co{Ms). Since A G Ac 
and Ac is open, A + eoC G Ac (here, e is of |y [-dimension) for some small enough eo > independent 
of S and v. Then z^ — [A + Coejs = {d — eo)e or [A + eoe]s = v + [eg — d)e. Since X + ege G Ac, 
A + eoe r^. Hence, eo - d < or d > e°. ■ 

The constant Cq will serve as a bound for the rate of the Lyapunov drift in the performance analysis. 

C. Performance Guarantee of LQF in Ac Region 

Therorem 12: If an arrival rate vector A satisfies A G Ac, then, the network is stable under the LQF 
policy. 

Sketch of Proof: Again, consider the fluid limit of the queue process and apply a similar argument 
as in the proof of Theorem [3l Let S* C y be the set of nodes whose queues are the longest at time t and 
will remain the longest for the next infinitesunally small time interval. Let us be the service rate vector 
for the nodes in 5 at time t. Under LQF, vs G Co{Ms) and us — [X]s = ee for some e. Since A G Ac, 
by Lemma \TT\ we have e > eo for some Cq > independent of S and v. Hence, at any time instance, 
each of the longest queues decreases at a positive rate no less than eo- This is sufficient to conclude that 
the original queueing process is a positive recurrent Markov process, which means the queues are stable. 
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Lemma 13: 0,° C Ac. 

Proof: Suppose Q° ^ Ac*. Then, there exists a vector A G 0° and A Ac. Hence, [X]s = u + de 
for some non-empty 5 C y, zy e Co{Ms) and d > 0. Since A G il°, A + ee G $7 for some small enough 
e > 0. From [A + ee]s = u + {d + e)e and d + e > d > 0, we have [A + ee]s > v- Hence, A + ee 0, 
leading to a contradiction. ■ 

Consider Example 2 in Section |TV-A| With the linear programming tools introduced in Section |Vl one 
can check that A^ G Ac but A^ Ac. This explains why A^ can be stabilized by LQF while A^ cannot, 
even though A^ > A^. 

D. Rank Condition and A/? Region 

For the same average arrival rate vector, whether the i.i.d. arrival processes have zero or non-zero 
variances leads to significantly different stability behavior (in the former case, the arrival processes are 
deterministic with constant rates). This issue has been discussed in jTH where the authors develop a 
queue separation result related to a rank condition about the matrices of the maximal independent sets. 
We next generalize the rank condition. Then, we extend Ac to a larger stability region Ar. We will 
show A/j can be stabilized under LQF when the arrival processes all have non-zero variances. 

Definition 7: Let 5 C y be a non-empty set. We call the matrix {Ms, e) the extended schedule matrix 
for S (or graph Gs)- Let R{Ms,e) denote the rank of the extended schedule matrix, i.e., the number 
of linearly independent columns in the matrix {Ms,e). We say S (or graph Gs) has a high rank if 
R{Ms,e) = \S\. Otherwise, we say 5 (or Gs) has a low rank. 

Suppose 5" C F is the set of nodes with the longest queues at some time instance. When the arrival 
has non-zero variances, the queue separation result suggests (Lemma 1 and Lemma 3 of f\M ): If the 
rank R{Ms) < \S\ — 2, then, with probability 1, the queue sizes of 5 will not stay identical in the next 
infinitesimal time interval. We find that the condition R{Ms) < 15*1 — 2 can be relaxed to R{Ms,e) < 
\S\ — 1, i.e., the low rank condition in Definition |7] The queue separation lemma (Lemma 1 of |fT4l) 
uses the assumption R{Ms) < IS*] — 2 to obtain a vector v such that v'e = and u' Ms = 0. Such a 
vector u still exists when the low rank condition in Definition |7] is satisfied. Then, every subsequent step 
in the proof of the queue separation lemma still holds. The low rank condition is a generalization since 
R{Ms) < |5| - 2 implies R{Ms,e) < \S\ - 1. 

Roughly speaking, when the variances are non-zero, the randomness in the arrival processes pressures 
the queues in S to move around in an |5| -dimensional space. This means that the |5| queues cannot be 
simultaneously the longest queues for a sustained period of time (in which case, the queue trajectory 
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moves along a line), unless the service can fully compensate the pressure from the arrival processes. 
But, full compensation is not possible in the low-rank case since the service rate vector lives in a lower- 
dimensional space. What will happen is that some subset of the queues in S with a high rank will dominate 
the rest. This is known as queue separation. The implication is that, in the case of non-zero variances, 
there is no need to consider the low-rank subsets of V when evaluating the performance degradation of 
LQF. The discussion motivates the following definition of A/j. 
Definition 8: The A/j region is defined by 

SCV,S with high rank 

In words, a vector A is outside Aji if and only if A G r5 for some node set S that has a high rank. 

By comparing the definitions of Ac and A/j, we have the following lemma. 

Lemma 14: Ac ^ Ar. 

In addition to the i.i.d and mutually independent assumptions, the following assumption on the arrival 
processes is needed for technical reasons (see |[T4l for their relevance). 

Al: (The large deviation bound on the anival processes) Let Ai{n) be the cumulative arrivals at queue 
I (at node I £ V) up to time n, and let A; be the average arrival rate at queue /. For each e > 0, 

pn^±J. - Xi\ > e) < /3exp(-n7(e)) for all n > 1, for some 7(e) > and ;9 > 0. 
n 

Therorem 15: Assume the condition in Al holds and assume the variance of the i.i.d arrival process 

to each node is non-zero but finite. If an arrival rate vector A satisfies A € Ar, then, the network is 

stable under the LQF policy. 

Sketch of Proof: Again, consider the fluid limit of the queue process and apply a similar argument 

as in the proof of Theorem |3] Let S C V he the set of nodes whose queues are the longest at time t 

and will remain the longest for the next infinitesimally small time interval. By replicating most of the 

arguments in the queue separation lemmas (Lemma 1 and Lemma 3 in HH), it can be shown that S 

must have a high ranly. Otherwise, the queue sizes of the nodes in S will be separated and they cannot 

all remain the longest. Hence, we can apply the same argument as that in Theorem [121 but only to the 

high-rank node sets. ■ 

■^The only change is to Lemma 3 in 1141 . Instead of saying for any low-rank set, there must be a subset that satisfies local 
pooling, we say for any low-rank set, there must be a subset that is of high rank. This is so because a set with a single node is 
of high rank. The modification is needed in the proof of Lemma 3 in 1141 . The statement of Lemma 3 also needs to be modified 
accordingly. 
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Some graph examples are given in Fig. |4l regarding their set cj-local pooUng factors and ranks. Note 
that, the shaded region includes those subsets S which either satisfy a^ = 1, i.e., set local pooling 
(SLoP), or have low rank. Those subsets need not to be considered for the performance of LQF in case 
of non-zero variances. 



Low Rank 


High Rank 


6-Cycle 


Petersen 


6-Cycle connected to 8-Cycle 


Mohius-Kantor 


SLoP 


Big-Cycle(\V\>=8) 
~~~~-^. Q^(Hypercube Graph) 


^--"^ 


/ Trees(\ V\ >2) 


Cliques Ny 


( 4-Cycle 


5-Cycle \ 


V^__ 


7 -Cycle J 





Fig. 4. Graph examples and classification by the set cr-local pooling factor and rank condition. For graphs Gs = (S, L) inside 
the oval, ag = 1; outside the oval, erg < 1. The graph labeled '6-cycle connected to 8-cycle' is shown in Fig. |6] 



E. Further Properties of Regions Ac and A/j 

It has been demonstrated that A^ CI Ajj. We now continue to study the properties of the two regions 
and their relationship. 

Therorem 16: The closures of Ac and Ar are the same, i.e. Ac = Ar. 

Proof: Since Ac ^ A/j, we have Ac ^ Ar. We will next show Ar C Ac- Since Ar = 



Ac [Ji^R - ^c) and Ar = Ac\J {Ar - Ac), we only need to show {Ar - Ac) C Ac- 
Given any vector A G Ar — Ac, by comparing Definition [6] and [8j we have 

^^ n r^^n( u r.). 



scv 

S with high rank 



SCV 
S with low rank 



By Lemma [TOl Ts is a closed set and F^^ is open. Hence, Ar is open. Therefore, there exists 6 > 
such that 7 S Ar whenever 7 > and the distance between the two vectors d{X, 7) < 6. 

■^^6 and A = A + ee. Then, the distance between A and A is d(X,\) = \ee\ < 
--6\e 



Let < e < „ , 

ij. Then, A G Ar and A > ee. 



2^\V\'''' 2 

Now, let (5(A) = {S\S <^ V,S with low rank, A G F^}. We will next construct a sequence of low-rank 
node sets. Si, for i = 1,2, •••. Since each of them has a low rank, there exists an |5i| -dimensional 
vector 5* / with ||5r*|| = 1 such that ((7*)'e = and {g'^)'Ms, = 0. We then extend each 5* to a 
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1 1/ 1 -dimensional vector by setting the values of the new components to be zero. With a little abuse of 
notation, we call this |y|-dimensional vector g^ as well. 

We now construct the sequence of Si. If Q{X) / 0, pick any subset ^i € Q(A). Let A^ = A + l/2eg^. 
Next, if Q{\^) / 0, pick any ^2 G Q{\^) and let A^ = A^ + l/2^eg'^. In step j, if Q{\^~^) / 0, we will 
pick any S^ G Q{\^~'^) and let A-' = A-'"^ + l/2^eg^. This procedure will go on until Q{X^) becomes 
empty for some j. We can check that the ith component of A-' is A-'(i) = (A + l/2eg^ + 1/2'^eg'^ + • • • + 
1/2^ eg^){i) > e — l/2e — l/2^e — • • • — l/2-'e > 0. This ensures that A-' is always a non-negative vector 
for all j. 

Now, we will show that there exists an integer Jv > such that Q{\^) becomes empty for the first 
time (hence, the sequence of Sj ends at Sk~i, or contains no elements ii K = 0). For convenience, let 
AO = A. 

We will show that Sj Q{\^) for k > j, where Sj / 0. Suppose Sj G Qi^'') for some k > j. 
Then, A'^ G Tg^, which implies that [A'^]^^ = d^e + u^ for some d^ > and v^ G Co{Ms-). From the 
construction procedure, we know that Sj G Q{\^~^), which implies that [\^~^]s, = d^e + v^ for some 
(i2 > and v'^ G Co{Ms^). Since 

A*^ = A^-^ + \l2^eg^ + \l2^^^eg3^^ + • • • + 1/2^/, 



we have 



Then, 



\\%^ = d'e + v' + [1/2^ eg^ + l/2^+'eg^+' + ■■■ + l/2'^eg%^ 



{g^)'[\%, ={9')'{d^e + v^ + [1/2^ eg^ + l/2^+\g^^^ + ■■■ + l/2hg%^) 

= l/2U\W\\^ + l/2^+^e{g^y[g^+%^ + • • • + l/2\{g^y[g%^ 

>l/2Je - 1/2^+^6 l/2''e > 0. 

However, since [A^^'J^^ = d^e + v^, we have {g^)'[\'^]s- = {g^y{d^e + i'^) = 0, leading to a contradiction. 
Hence, Sj Q{X'') for k > j. 

In summary, each non-empty Sj in the constructed sequence is in Q{X^~^) but not in QiX'') for k > j. 
Hence, each Sj is distinct. Since there is a finite number of non-empty node sets S C y, there exists an 
integer A' > such that Q{X^) becomes empty for the first time. 

Then, A^^ Fg for any node set S with a low rank. Hence, X^ G Hscyswithiowrank-'^s- "^^^ distance 
between A and A^*^' is d{X, X^) < e(l/2 + 1/2^ + ... + 1/2^) < e. Then, the distance between A and A^ 
is d{\X^) < d(A,A) + d(A,A^) < e^/\V\+e < 5. Hence, A^ G A/?, It follows A^ G nscyTs = A^. 
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Since e can be chosen arbitrarily small, A is a limit point of Ac. Thus, A G Ac, implying (A/j — Ac) ^ 



Ac- Hence, {Ar - Ac) C Ac- ■ 

The following is an intermediary lemma. 

Lemma 17: If a non-empty set 5 C F satisfies a'^ = 1, then Fg P| A° = 0. 

Proof: Suppose there exists a vector A G Tg P| A°. By Definition \5\ [X\s = i^ + de for some d > 
and u G Co{Ms). Since A G A°, by Lemma [T] [Ajs + ee < /U for some fi G Co{Ms) and a small enough 
e > 0. Hence, fi> u + {d + e)e > z/(l + d + e). Thus, u^ < 1 and we arrive at a contradiction. ■ 

Lemma 18: If every high-rank node set 5 C y satisfies a^ = 1, then, Ac = A/j = A and Ar = A°. 
Proof: According to Definition [H we have Ar = Clscv s withhi'yhmnk^s- ^^^ ^^y high-rank node set 
S, since a^ = 1, we have T^ fl A° = ^ by Lemma [HI which impUes T^ fl ^° = ^°- Hence, Ar f]A° = 
r\scv,s with high rank T^ fl A° = ^"- Combining this with Theorem [HI we get A^ = A^ = A^ = A. Also, 
the fact that Ar f] A° = A° implies A° C Ar. Since Ar is an open set in A and A° is the largest open 
set in A, it must be that Ar = A°. ■ 

Remark. From Lemma [T8l we know that when all the subsets S* C y satisfy either set local pooling (i.e., 
ag = 1) or the rank of S is low, then Ar = A°. That is, the entire A" is achievable by LQF, assuming 
the arrival processes have non-zero variances. This is the same statement as Theorem 1 of |[T4l . Thus, 
the newly developed theory here is able to reproduce the main result of |[T4ll . 

Lemma 19: Ac = A if and only if Ac is convex. Similarly, Ar = A if and only if Ar is convex. 
Proof: It is obvious that Ac = A implies Ac is convex. We will next show the converse. Since 
Ac ^ A and A is a closed set, we have Ac ^ A. Because 0,° C Ac, we have U C Ac- Since U 
contains all the extreme points of A (Corollary |7]), Ac also contains all of them. Since A is the convex 
combination of all its extreme points and Ac is convex, we must have A C Ac. 

The second statement can be proved similarly. ■ 

V. Graph coloring and LQF scheduling 

The scheduUng problem in this paper is deeply connected with graph coloring and its related problems. 
In this section, we will introduce fractional coloring, and more generally, aspects of the fractional graph 
theory that can provide useful tools for studying the stability regions discussed in the previous sections. 

A. Fractional Coloring and Capacity Region 

The chromatic number of a graph G, denoted by x{G), is the minimum number of colors needed to 
paint the nodes so that the connected nodes do not share the same color. When we relax the integrality 



20 

constraints of the chromatic number problem and introduce a parameter A S m1^ , we obtain the following 
linear programming (LP) problem. 

Definition 9: Given a graph G = {V, E) and A G m1^ , the weighted fractional coloring problem with 
the weight vector A is: 

Xf{G,\) = min e'a, subject to Mya > A,a > 0. (9) 

The optimal value of the above problem, Xf{G, A), is called the weighted fractional chromatic number, 
which is known to be related to the capacity region as follows (see ||2TI ): 

A = {AgM^^I I X/(G,A)<1}. (10) 

Based on Q, Xf{G,X) can be interpreted as the fastest way of serving queued data when the queue 
sizes are proportional to the weights A. Based on (ITOl ). XfiG,X) can be interpreted as the 'traffic load' 
to the network. 

The relevance and usefulness of this problem to the study of wireless scheduling have been amply 
demonstrated in lIlTI . The characterization of the capacity region by ( fTOl ) suggests that the fractional 
chromatic number can serve as an oracle for judging whether an arrival rate vector is in the capacity 
region or not. With this observation and with known complexity results about the fractional coloring 
problem, the authors of llTl have derived results about the inherent complexity of the wireless scheduling 
problem. 

B. Weighted Fractional Matching Number and il. Region 

We next discuss the problem of finding the weighted fractional matching number of a graph |[T9l . This 
problem can help to decide whether a vector is in Q,. 

Definition 10: Given a graph G = {V,E) and A G M!^ , the weighted fractional matching number 
problem with the weight vector A is: 

(j)j{G, A) = max e'/3, subject to My/3 < A,/3 > 0. (11) 

The above problem is the Lagrangian dual of the weighted fractional transversal number problem, 
which is the hypergraph dual problem of the weighted fractional coloring problem |[T9l . Here, the ith 
component of /3 can be interpreted as the amount of time for which the ith maximal schedule is used. 
The weighted fractional matching number, (/)j(G, A), can be interpreted as the slowest way of serving 
the queued data (in the amount A) using only the maximal schedules, subject to the additional constraint 
that a schedule should not be selected if it activates a link associated with an empty queue. 
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Lemma 20: The $7 region satisfies the following: 

= {A G M^^l I 0/(G5, [A]5) < 1,V5 C y,5 / 0}. 

Proof: Consider any vector A G J7 and an arbitrary non-empty node set 5" C F. Suppose (l)f{Gs, [^]s) 
> 1. Then, by Definition [TOl we have [X]s > ku for some z/ E Co(M5') and A: > 1. Let Z be the largest 
subset in 5 such that [i^]z = 0. Note that Z ^ S. Then, the vector [i^js^z G Co{Ms^z) and [i^]5-z > 0. 
Hence, [AJ^-z > fc[i^]s'-z > Ms-z, which implies that A G Hs-z- According to Definition [3] and IH 

Conversely, suppose a vector A satisfies (pj^Gs, [X]s) < 1 for every non-empty S CV. Then, [A]^ ;2i z^ 
for any u G Co{Ms). Otherwise, there would exist k > 1 such that [X]s > ku for some z^ G (70(71^5), 
which implies that cl)f{Gs, [X]s) > 1- Thus, A G f]. ■ 

C. Hypergraph Duality and Set a-local Pooling 

We next relate the ratio of Xf{G, A) and (pf{G, A) to the set a-local pooling factor. First, we have the 
following lemma. 

Lemma 21: For any /c > 0, Xf{G, kX) = kxf{G, A) and (pfiG, kX) = kcl)f{G, A). 

Proof: The case of /c = is trivial. We only focus on the case of fc > 0. Suppose /3* is an optimal 
solution to the problem in (ITTI ) for finding (j)f{G,X). Then, k(3* is feasible to the problem for finding 
(t)f{G,kX). Since e'{kfi*) = ke'/3*, (pf{G,kX) > k(Pf{G,X). 

Conversely, suppose /3 is an optimal solution to the problem for finding (l)f{G, kX). Then, f3/k is feasible 
to the problem for finding (j)f{G,X). Since e'{(3/k) = efi/k, (l)f{G,X) > (j)f{G,kX)/k. Therefore, 
(j)f{G,kX) = k(l)f{G,X), for A; > 0. A similar argument can be used to show Xf{G,kX) = kxf{G,X). 

m 

Therorem 22: Given a non-empty node set S C V, the set u-local pooling factor of S satisfies the 

following^ : 

* . X/(Gs,A) 

as = mm -7^-— — — . (12) 

x>o (l)f{Gs,X) 

^We take the convention a/0 = oo for any scalar a > 0. Note that, if any component of A is equal to zero, then (f>f{Gs, A) — 0. 
As a result, the optimal solution A* to ( I12t must satisfy A* > 0. 
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Proof: Suppose (a^, ^^, Vg) is an optimal solution to the problem in (|5]l. Then, we choose A = Vg. 
Since o"|./i^ > z^^, we have 

Xf{Gs,\) = Xf{Gs,y%) 

< XfiGs,cr*sHs) = (J*sXfiGs,fJ'*s) < CTs- 

The last inequality above uses the fact Xf{Gs, fJ-g) < 1, which follows from ( fTOl ) (since Hg G As). 

Because Ug G Co{Ms), there exists a non-negative vector /3 such that M^/S = i/^ and e'/3 = 1. Such /3 
is feasible to (fTTI ) for finding (f)f{Gs, A). Hence, (j)f{Gs, A) > 1. Therefore, minA>o Xf{Gs, ^)/4'f{Gs, ^) < 

Next, suppose A* is an optimal solution for the problem minA>o X/l^s, A)/(/>/(G5, A). Suppose 
a* and /3* are optimal solutions for the problems of finding Xf{Gsi^*) and (f)f{Gs, X*), respec- 
tively. Then, we have m.mx>oXf{Gs, X)/<Pf{Gs, X) = E^»i/Ei(^t■ Now, let fi = Msa*/Y.i0^t 
and V = MsP*/ J2i P*- Then, fi,u e Co{Ms) and 

Xf{Gs,X*) _ T.i<Msa* _ Msa* 

By the feasibility of a* and j3* to ^ and ([Tl]), respectively, M^a* > A* > Ms (3*. Hence, 

Msa* ^ X* ^ MsP* _ 

Thus, XfiGs, X*)/(t)f{Gs, A*) is feasible to the problem in ^. Therefore, minA>o Xf{Gs, X)/cl)f{Gs, A) > 
a*s. m 

The theorem above shows that the set cr-local pooling factor is the same as the minimum of the 
hypergraph duality ratios over all different weights. 

D. Weighted Fractional Domination Number and Ac Region 

Definition 11: Given a graph G = {V, E) and A G W_^_ , the weighted fractional domination number 
problem with the weight vector A is: 

Tf{G, X) = maxd, subject to de + v = X,u £ Co{Mv)- 

For convenience, let tj(G, A) = — oo when the problem is infeasible. 
By Definition |5l [6] and [TT] we have the following lemma. 
Lemma 23: The following relations hold: 

Ts = {Xe M^^l I Tf{Gs, [X]s) >0, for S<ZV,S^ 0}, 

Ac = {A G M^^^l I TfiGs, [X]s) < 0,V5 CV,S^ 0}. 
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VI. Experimental Examples 

In this section, we show some simulation results. The main purpose is to confirm some of the less 
intuitive theoretical results. We first show the performance of LQF on the six-cycle graph, denoted by 
Cq, for aiTival rate vectors in different sections of the capacity region. For Cq, LQF can achieve the 
entire interior of the capacity region for arrivals satisfying assumption Al and with non-zero variances. 
On the other hand, for constant arrivals, experiments have shown that some rate vectors in the interior 
of the capacity region are not achievable by LQF. 

In the experiments with constant arrivals, we use a load parameter to scale the arrival rate vectors. 
Each experiment runs for 10^ iterations with initial queue sizes of 10^. The following arrival rate vectors 
are used for the results in Fig. [S) 

,1 1 1 1 1 1 

A = ( e, e, e, e, e, e) 

^2 2 2 2 2 2 

,,11 1 11 1 ^, 

^ =(2-^'2-'^'2-^^'2-^'2-'^'2-^^) 

o ,1 1 1 1 1 1 X, 

^ =(2-^'2-'^'2-2^'2-2^'2-'^'2-2^)' 
where e = 10~^. Note that 0.7A^ < 0.95A^ < A'^. However, judging by the queue sizes in Fig. |5l the 
arrival rate vectors 0.7A^ and 0.95A^ seem to be not stabilizable, whereas A'^ seems to be stabilizable. 
The theory allows this counter-intuitive phenomenon. Readers can verify that 0.7A^,0.95A^ ^ Ac while 
A=^ G Ac. 
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Fig. 5. Constant arrivals in Ce 



In Section |lVl we generalize the definitions of high or low-rank graphs. In Fig. [6j we provide an 
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interference graph that is not set local pooling (crj < 1) and is of high rank according to the original 
definition in |[T4ll . However, in the new definition, the graph is of low rank. We ran simulations with 
initial queue sizes of 10^ using the Bernoulli arrivals with an identical arrival rate of 0.499. Fig. |7] shows 
the evolution of the average queue size for nodes 1-8 and 9-13 over 10^ iterations. The queues for nodes 
1-8 appear to be unstable and the queues for nodes 9-13 appear to be stable. Our refinement of the rank 
condition rules out the possibility that the queues of all nodes are simultaneously the longest and remain 
longest, whereas the previous rank condition does not rule that out. 




Fig. 6. An interference grapli with Ca connected to Cg. For this graph, the ranks are R{Mv) ~ R{Mv, e) = 12. 
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Fig. 7. Average queue sizes for nodes 1-8 on Cs (labelled as Vji^gj) and for the nodes 9-13 on Ce (^[9,12]). 



VII. Conclusion 

In this paper, we investigate the performance guarantee of the LQF scheduling policy in wireless 
networks. The objective is to discover new stability regions of LQF that are larger than those previously 
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known, and to improve our knowledge about the largest possible stability region of LQF. We show that 
it is necessary to go beyond the existing framework of linear reduction of the capacity region, and move 
to a non-linear framework. 

We introduce the concepts of strictly dominating vectors and uniformly dominating vectors; the former 
leads to the new stability region of LQF, Q,, and the latter leads the stability regions Ac and A/j. We 
show that O contains S*(G)A, which is the stability region given in ifTSl . We also show $1° C Ac ^ A/j. 
Hence, the new stability regions all capture the performance of LQF better. Contrary to the previously- 
known regions of stability, the closures of these new stability regions contain all the extreme points of the 
capacity region A, but they are not convex in general. The only case where they are convex is when they 
are equal to the capacity region itself, which occurs only for selected interference graphs. The general 
lack of convexity is not a defect of the theory. We show that, when LQF cannot achieve the full capacity 
region, the largest achievable region cannot be convex. 

The study reveals a counter-intuitive situation where increasing the arrival rates helps LQF to stabilize 
the network. It turns out, in this case, the original rate vector is outside Ac, and after the rate increase, 
the new rate vector is inside Ac. We also generaUze the rank condition studied in |[T4l . and with this 
generalization, refine the stability results for non-deterministic arrivals. We can show that if a set of nodes 
satisfies the new low-rank condition, the queue sizes of these nodes will be separated. Based on this result, 
we can enlarge Ac to Aji, which is achievable by LQF under non-deterministic arrivals. Interestingly, 
we show that the closures of Ac and Ar are the same. Finally, we introduce several linear programming 
problems encountered in the fractional graph theory, which can provide tools for studying the newly 
developed stability regions. We show that a ratio between the weighted fractional coloring number and 
the weighted fractional matching number is related to the set cr-local pooling factor introduced in lITSl . 

VIII. Appendix 

Proof of Lemma \5l[ Suppose u is an independent set of G, represented by a 0-1 vector. Clearly, 
1/ G A. Let us write u = av^ + (1 — cl)v'^ for some < a < 1, and u^, v^ G A. Note that /i < e for every 
[I € Co{},Iy\ Thus, < z^^ < e and ^ < v^ < e. For any index i, if vi = 0, we must have v\ = and 
v1 = 0. Similarly, if i/j = 1, we have uf = uf = 1. Therefore, v = v^ = u'^ and v is an extreme point 
of A. 

Conversely, take any extreme point v of A. Then, v < fi ioi some fi G Co{Mv)- For an index i, if 
Ui > 0, we claim that Vi = fii. Otherwise, we let t = fii — Vi > 0. Then, we can create P and D such 
that Uj = Vj = Vj for j ^ i. We can find an e > such that Vi = Vi — et > 0, and we let Vi = Hi > 0. 
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Since < P < P < /^, we have P, z> G A. It is easy to see v = YirJ^ + tT^^' which implies v is not an 
extreme point. Hence, either Vi = or Vi = fii, for all i. 

We now show that Vi = or Vi = fii = 1 for all i. Write fi as fi = "^'j^i ajfi^ , where each fi^ G My, 
each Qj > 0, and Xlfci '^j — 1- Let f^ = if f j = 0; f^ = /u;^ otherwise, for all < j < k. Because 
fi^ G My, we have t'^ = or vf = 1. It is easy to check that u = X]7=i o,jiy^. Since z^-' < fi^, we have 
z^-' G A for each j. Since u is an extreme point of A, u^ = u^ = ■ ■ ■ = u'^. Thus, Ui = 1 or i>i = for 
all i 

It is easy to see that any 0-1 vector in Co{Mv) must be a feasible schedule, i.e., an independent set 
of G. Since fi G Co{Mv), the set of nodes, S, for which fii = 1 forms an independent set. Let 5' be 
the set of nodes for which z/j = 1. We have S' C S. Therefore, u corresponds to an independent set. ■ 
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